Charge order, dynamics, and magneto-structural transition in multiferroic LuFe204 
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We investigated the series of temperature and field-driven transitions in LuFe2 04 by optical and 
Mossbauer spectroscopies, magnetization, and x-ray scattering in order to understand the interplay 
between charge, structure, and magnetism in this multiferroic material. We demonstrate that 
charge fluctuation has an onset well below the charge ordering transition, supporting the "order 
by fluctuation" mechanism for the development of charge order superstructure. Bragg splitting and 
large magneto optical contrast suggest a low temperature monoclinic distortion that can be driven 
by both temperature and magnetic field. 
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Complex oxides take advantage of the unique and flexi- 
ble properties of transition metal centers to govern bond- 
ing and local structure. Further, the delicate interplay 
between charge, structure, and magnetism yields impor- 
tant consequences for functionality and cross-coupling. 
Iron-based materials such as multiferroic BiFeOa and 
LuFe204 P, bistable photomagnetic systems such as 
Prussian blues and related derivatives 3 , dilute magnetic 
semiconductors Q, and the new family LaFeAsOi-^^Fa, 
of superconductors [Hi have attracted recent attention. 
In this Letter, we focus on LuFe2 04, a frustrated sys- 
tem with a series of phase transitions that give rise to 
electronically-driven multiferroicity Q . 

LuFe204 has a layered structure with Fe-containing 
double layers of triangular connectivity. Three- 
dimensional Fe^+/Fe'^+ charge order (CO) occurs below 
320 K (Tco)- This is followed by ferrimagnetic order be- 
low 240 K (Tn) 0,0, H- An additional low temperature 
magnetic phase transition has recently been reported at 
175 K (Tlt) S|- The CO has a so-called \/3xV^ super- 
structure 0, U Because of the mixed valcnt iron centers 
and frustrated triangular lattice, the Fc^+ and Fc^+ pop- 
ulations are different within the double layer, an effect 
that renders an intrinsic polarization 0, @| . The charge 
ordering mechanism is thus central to understanding the 
unusual physical properties of this multiferroic. 

To elucidate the charge excitations and understand 
how they correlate with structure and magnetism, we 
measured the optical and Mossbauer spectra, magnetiza- 
tion, and x-ray scattering of LuFe204. We compare our 
comprehensive results to recent electronic structure cal- 
culations 01 and to spectral data on classical magnetite 
[lot . We demonstrate that strong Fe^"*" ^Fe^+ charge 
fluctuation persists even in charge ordered states charac- 
terized by superstructure reflections, and it persists down 



to Tlt below which Bragg splitting indicates that strong 
monoclinic distortions occur. These observations are con- 
sistent with the "order by fluctuation" mechanism , in 
which case -x/Sx a/S CO is preferable for entropy reasons 
and stabilized by the charge fluctuation in this geomet- 
rically frustrated system. As in magnetite, we analyze 
the results in terms of a polaron picture, extracting a 
large effective mass for the charge carriers. On the other 
hand, Fe^+ on-site excitations are sensitive to the mag- 
netic transition at Tlx and display a sizable magneto- 
optical effect. Combining our spectral, magnetic, and 
structural data, we generate an H-T phase diagram and 
show that the transition at Ti^t can also be driven by a 
magnetic field. These results demonstrate the important 
interplay between charge, structure, and magnetism. 

All experiments were conducted on floating-zone- 
grown LuFe204 single crystals from the same batch as 
those used in Refs. [1,01. Near-normal reflectance mea- 
surements were carried out on ab plane samples employ- 
ing a series of spectrometers covering a wide range of 
energy (30 meV - 6.5 eV), temperature (4 - 540 K) and 
magnetic field (0 - 33 T, 1| c) [l21- Optical conducti vity 
(Ji{E) was calculated by a Kramers-Kronig analysis [l3|. 
Variable temperature transmittance was done on a 25 ^m 
ab plane crystal, allowing direct calculation of absorption 
a{E). The Fe-57 Mossbauer spectra of 35 mg/cm^ of 
crushed crystals were recorded between 260 and 400 K 
on a constant acceleration spectrometer with a Rh matrix 
Co-57 source and calibrated at 295 K with a-Fe powder. 
The reported isomer shifts are relative to a-Fe at 295 
K. X-ray scattering was performed at undulator beam- 
line 4-ID-D at the Advanced Photon Source with 36 keV 
photons employing a cryostat with a split coil vertical 
field magnet up to 4 T. The sample was mounted with 
an angle of 45° between c and the field as a compromise 
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FIG. 1: (Color online) (a) Mossbauer spectra of LuFe204 and 
fits, see text; bottom: alternative best fit with a constraint of 
33% of hopping electrons (b) Arrhenius plot of the hopping 
frequency. Open symbols correspond to the spectra in (a), 
bars to hopping frequencis in [ij] . (c) Mossbauer spectral pa- 
rameters, from top to bottom: hopping frequency, quadrupole 
splitting, isomer shift, and full linewidth at half maximum. 



between cryomagnct angular restrictions and accessibil- 
ity of important regions of reciprocal space. Previous 
magnetization work shows that LuFe204 is rather insen- 
sitive to fields II ah below 7 T. Consequently the dominate 
effect of the field is due to the component || c. 

The Mossbauer spectra, (Fig. [TJa)) were fit with a 
Blume-Tjon model [15| for Fe^"*" and Fe"^"*" relaxation, 
similar to [l3|. individual fits of the spectra reveal (i) 
two Arrhenius processes as indicated by the temperature 
dependence of the hopping frequency and break at Tqo 
(Fig. mb)), (ii) the difference between Fe^+ and Fe^+ 
isomer shifts, Ai5, is constant below Tqo and cannot be 
resolved above, and (iii) a constant linewidth and a grad- 
ual ~ \JTco — T broadening above and below Tco, re- 
spectively. In order to overcome correlation effects (see 
below) and to reduce the large number of fit parameters, 
a simultaneous parametric fit of all spectra was then car- 
ried out with the constraints (i), (ii) constant A(5 for all 
T, and (iii). The fits in Fig. [DJa) are the result of this 
simultaneous fit. The obtained spectral parameters and 
relaxation frequencies (Fig. [ljb,c)) are in agreement with 
14| . with the exception of the sharper T dependence of 
these parameters around Tqo i which we attribute to the 
preparation of the sample as single crystal. Above Too, 
we obtained an activation energy of 0.16(2) eV. This is 
somewhat smaller than the 0.26 eV e nerg y obtained by 
electrical conductivity measurements [1^. Below Tco, 
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FIG. 2: (Color online) Optical properties of LuFe204. (a) 
(Ji vs energy E at 540 and 60 K calculated from a Kramers- 
Kronig analysis of reflectance. Brackets indicate assignments, 
the hierarchy determined from Ref. 0] . (b) a ys E calculated 
from transmittance, showing the charge transfer edge starting 
from ~0.5 eV. Inset: V aE with linear fit (dotted line). (c,d) 
ne//(r) calculated by sum rules for Fe'^^ — >Fe^^ charge trans- 
fer (c) and Fe^^ on-site (d) excitations, (e) Eg{T) calculated 
from an absorption edge fit assuming an indirect gap. 



we find the activation energy to be 0.36(9) eV. The large 
error bar is due to (i) small hopping frequencies, close to 
the detection limit, and (ii) correlations between the hop- 
ping frequency and linewidth. The Fe^^ line around 1.4 
mm/s is increasingly broadened below Tco due to micro- 
scopic lattice distortions in the charge ordered state (Fig. 
IHa)) [11. The fit at the bottom of Fig. [B^a) is the best 
fit with a constraint of 33% of Fe^-^+ with a hopping fre- 
quency of 1 GHz. The poor fit quality indicates that the 
260 K hopping frequency is smaller than 1 GHz, which 
invalidates the presence of 33% of Fe^-^"*" below Tco sug- 
gested in Ref. 1^ from modeling electron diffraction, a 
technique with a resolution better than 1 GHz. 

To further study the charge fluctuation (Fe^+ ^Fe'^+ 
charge transfer) observed in Mossbauer spectroscopy, we 
employed optical spectroscopy, a technique in which hop- 
ping is driven by comparatively high frequency photons 
(^ 10^^ Hz, in contrast to the natural hopping rate of a 
few MHz, as observed by Mossbauer spectroscopy). Fig- 
ure m^a) displays the optical conductivity of LuFe2 04. 
We assign the observed excitations based upon recent 
first principle calculations 0. The lowest allowed elec- 
tronic features are minority channel Fe^'*' — s-Fe^"*" charge 
transfer excitations. At slightly higher energy follow Fc^+ 
on-site excitations. Minority channel O p-^Fe d charge 
transfer and overlapping majority channel O p ^ Lu s 
state excitations are observed at higher energy (>3 eV). 
It is difhcult to resolve all the excitations because they 
are broad and overlap significantly. Nevertheless, the 1.1 
eV peak in the near infrared range (Fig. [2Ia)) can be as- 
sociated with Fe^+ ^Fe'^+ charge transfer. The optical 
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gap Eg is determined by the absorption edge of this band. 
A close-up view of the tail is shown in Fig. [UJb). The 
linear fit of \/aE vs. E above 0.5 eV indicates that the 
gap is indirect in agreement with Ref. [l^. Due to 
the 0.3 eV feature (possibly a spin-forbidden Fe'^"'" on-site 
excitation [111), the typical "double slope character" is 
not observed and the associated coupling phonon energy 
for the indirect gap excitation process can not be deter- 
mined. The optical gap is sensitive to Tco (Fig. E^e)), 
although LuFe204 is a semiconductor (non-metallic) over 
the full temperature range of our investigation. 

To quantify the strength of the various excitations, we 
calculated the effective number of electrons Ueff from 
the optical conductivity (Ti(a;) using the partial sum rule: 

■n-eff = J CTi(w)/eodw ^ ^ttuj^ , where ujp = ^ 

is the plasma frequency, e and m are the charge and mass 
of an electron, eo is the vacuum dielectric constant, Vq 
is the unit cell volume, and uji and uj2 are the frequency 
limits of integration. The absolute level of Ueff depends 
on the integration range, but the temperature trends are 
not sensitive to this choice. For instance, to investigate 
changes in the Fe charge transfer band, we evaluated the 
partial sum rule from 0.6 - 1 eV, as indicated in Fig.[5Ka). 
In this case, rie/ / represents effective number of electrons 
that are able to overcome the energy barrier to hop from 
Fe^"'' to Fe"^+ sites. This number increases over a broad 
temperature range through Tco, a-s shown in Fig.[2ljc). 



In optical processes, electrons hopping from Fe + to 
Fe'^^ are better described as small polarons, which cor- 
respond to combined electronic and vibrational excita- 
tions that arise when the lattice is too slow to relax 
[lit . Important signatures include (i) a large effective 
mass and (ii) optical excitation energies that are much 
larger than the low frequency activation energy. The ef- 
fective mass of the charge carriers can be estimated us- 
ing n^ff = where to* is the effective mass and 
= 3 (the number of Fe^+ site per unit cell) . Consid- 
ering that the 0.6 - 1 eV integration is only over half of 
the excitation, we get ^ w 40, which is large, but typ- 
ical for polarons (e.g. ^ « 100 in Fe304 0). With 
the polaron picture and the simple model of an electron 
jumping between two sites [11], we can estimate the 60 
K activation energy from the optical activation Hloq as 
Ea = fttJo/4=l.l(l)eV/4=0.28(3)eV, in excellent agree- 
ment with value reported in the study of low frequency 
dielectric dispersion and DC electric conductivity rang- 
ing from 0.25 - 0.29 eV [2, 14| and compatible with that 
for spontaneous electron hopping obtained between 260 
and 320 K from Mossbauer spectroscopy. 

The T dependence of Ueff corresponding to Fe^+ to 
Fe'^+ charge transfer confirms charge fluctuation below 
Tco (Fig. 2(c)). Here, Uefj begins to increase well 
below Tco (evident also in electron hopping trends via 
Mossbauer spectroscopy) and continues to change above 
this temperature. This result is consistent with the 



16 
14 
12 
10 



1 

o 



o 1.4 
ri 

1.0 







4.2 K 
550 nm 







H\ 




Optics: 

-□- R(H ), H t 
Magnetization: 

M(H ),Ht 
-^r— M{H ),Hl 
—*-M(T). Ti 
X-ray scattering: 
o{336) not split 
■ •(336) split 
<) (336) split/not split for Hm 







50 



100 



150 r(K) 



FIG. 3: (Color online) Low-temperature phase transition 
from magnetization data (AT -4, see [^), field-dependence of 
optical reflectance contrast R{H)/R{H = 0) at A = 550 nm 
(squares, see left inset), and presence (indicated by •) or ab- 
sence (o) of a splitting in the total diffraction angle of the 
structural Bragg reflection (336) from synchrotron x-ray scat- 
tering (right insets). The maximum region of the structurally 
distorted phase is shaded (see text). 



presence of 3D anti-phase domain boundary modes 21 1 . 
We attribute the experimental observation of significant 
charge fluctuations through Tco (even where diffraction 
shows that it is ordered) to relevance of the "order by 
fluctuation" mechanism in which fluctuations are needed 
to stabilize the \/SxV3 CO in the frustrated system [ll[. 
Interestingly, the charge fluctuation onset is at Tlx (Fig. 
2(c), (e)), suggesting that the low temperature phase 
transition quenches the charge fluctuation. Similar to 
the Verwcy transition in FC3O4 [23, the T dependence 



of rieff shows an anomaly near Tco consistent with the 
lowering of the activation energy above Tco [3l- How- 
ever, the jump occurs above Tco- 

The strength of the Fe^+ on-site crystal field excitation 
is quantified by the partial sum rule in energy range 2 - 
3 eV, in accord with first principle calculations Al- 
though the data are more scattered than that discussed 
above because of the background from nearby excita- 
tions, two finding are immediately clear. First, the Fe^+ 
on-site excitation is rather insensitive to charge and spin 
ordering transitions. Second, it displays a clear anomaly 
at Tlt (Pig-Wic)), which recent magnetization and neu- 
tron diffraction studies Q identified as an additional first- 
order transition. Further magnetization measurements in 

||c up to 7 T indicate that the transition temperature is 
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strongly suppressed by H (Fig. [3]). Hysteresis of Hlt{T) 
widens upon cooling, and below 50 K the high-i/, higli- 
T phase remains frozen in even after decreasing 7J to 0, 
reminiscent of kinetic arrest of first-order transitions as 
studied, e.g., in doped CeFe2 psj . 

To check for possible Q structural components of the 
transition, we closely examined structural Bragg reflec- 
tions, such as (336). CO below 320 K [§] lowers the 
crystal symmetry to monoclinic, which could lead to a 
monoclinic distortion {(3 ^ 90°), with a splitting in 29 
values of such reflections, with domain formation clearly 
observable also in single crystals [l^. Such a splitting 
could not be resolved in between 200 and 300 K (in any 
magnetic field) , indicating that any monoclinic distortion 
is small. Splitting becomes evident upon cooling below 
Tlt (lower right inset in Fig. [3]), consistent with a signif- 
icant (/3 « 90.5°) monoclinic distortion [2^. Application 
of a magnetic field removes the splitting (upper right in- 
set) as soon as the field component ||c reaches Hlt as 
determined by magnetization, and subsequent reduction 
causes a reappearance of the splitting at a field value 
again consistent with magnetization. These diffraction 
data thus suggest that Tlt and Hlt{T) corresponds to 
a strongly hysteretic magneto-structural transition. We 
propose that the monoclinic distortion removes geomet- 
ric frustration rendering charge fluctuation unnecessary. 
This scenario is in line with the observed fluctuation on- 
set at Tlt and the "order by fluctuation" mechanism. 

The anomaly at Tlt in the Fe^+ on-site excitations can 
also be explained within a structural distortion scenario. 
Consider an Fe center coordinated by five O ligands in a 
trigonal bipyramidal geometry {D^d symmetry) . A crys- 
tal field splits the Fe 3d levels into three groups [lH . The 
monoclinic distortion in the low T phase splits these lev- 
els further, shifting the on-site excitation energies and 
causing the discontinuity in rie// (Fig. EJd)). 

The upper-left inset in Fig. [3] displays the 4.2 K re- 
flectance of LuFe204 at 550 nm as a function of mag- 
netic field {H |j c). At this energy, the spectral response 
is probing field-induced changes in Fe^+ on-site excita- 
tions. Strikingly, the reflectance increases by ~40% for 
H > lAT. This is consistent with a straight-forward ex- 
trapolation of Hlt from the magnetization data (dashed 
line in Fig. [3]). That the reflectance maintains this high 
value even after subsequent complete removal of H is 
again consistent with the magnetization data (sec above) . 
Thus, at low T a rather large field is required to switch 
the crystal structure. Since the structural distortion is 
most likely induced by the CO, which lowers the crystal 
symmetry Q , the strong i7-dependence of this transition 
is a further example of the strong coupling between spin, 
structural, and charge degrees of freedom in LuFe204. 

In summary, optical and Mossbauer spectroscopies 
demonstrate that charge fluctuation in LUFC2O4 has an 
onset at Tlt, well below Tco, supporting the "order by 



fluctuation" mechanism for the -\/3 x V3 CO superstruc- 
ture. Fe^+ on-site crystal field excitations are sensitive 
to the magneto-structural transition at Tlt, which can 
be driven by both temperature and magnetic field (re- 
quiring 14 T at 4 K). Combining spectral, magnetic, and 
structural data, we generate a comprehensive H-T phase 
diagram. The large temperature range of the observed 
dynamical effects is a consequence of the strong coupling 
between charge, structure, and magnetism. 
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